Nature of mechanical instabilities and their effect on kinetic friction 
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It has long been recognized that the key to understand kinetic friction force Fk is the analysis of 
microscopic instabilities that lead to sudden irreversible "pops" of certain degrees of freedom. In this 
Letter, the nature of such instabilities is characterized with an emphasis on boundary lubricants. It is 
shown that there are certain critical values of the parameters defining our model Hamiltonian, where 
the behavior of the instabilities changes qualitatively. Simultaneously, the functional dependence of 
Fk on the sliding velocity vq changes. The relevant parameters studied here are dimensionality of 
the interface, degree of commensurability, first higher harmonic in the lubricant wall potential, and 
temperature. Molecular dynamics simulations are carried out to test whether the predictions made 
on the basis of the simple model also hold in less idealized circumstances. 
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When a solid slider is moved laterally with respect to 
a substrate, the kinetic friction Fk is usually almost in- 
dependent of the sliding velocity vq jjj] with leading cor- 
rections in the order of In vq || . This so-called Coulomb 
friction differs from drag or Stokes friction that states a 
linear relation between Fk and vq and that can be under- 
stood from equilibrium statistical mechanics as a special 
case of the fluctuation dissipation theorem [||: For in- 
stance, the drag force experienced by a Brownian particle 
in solution arises from the many collisions between the 
Brownian particle and the solvent molecules. In general, 
the condition for a linear friction force to occur is the (an- 
harmonic) coupling of a central degree of freedom like the 
Brownian particle or a phonon to (infinitely) many other 
discrete degrees of freedom. In this Letter, the question 
will be addressed what the condition is for Coulomb fric- 
tion to occur. 

It has long been recognized that Coulomb friction must 
be related to instabilities that occur on a microscopic 
scale 0. When a slider is moved laterally with respect 
to the substrate, fast motion (pops) of certain degrees 
of freedom become unavoidable even if the sliders cen- 
ter of mass velocity is extremely small. The microscopic 
peak velocities in such pops are rather independent of Vq 
and consequently the energy dissipated via a Stokes-type 
mechanism also becomes almost independent of v$. The 
prototypical instability leading to Coulomb friction was 
suggested by Prandtl and Tomlinson H. In their one- 
dimensional model, a surface atom in the slider is coupled 
with a spring of stiffness fc to its ideal lattice site which 
moves at constant velocity vq- Interactions with the rigid 
substrate are modeled with a potential energy surface V 
that is periodic in the substrate's lattice constant plus 
some drag force linear in the atom's veloctity x. If k is 
sufficiently small, pops become unavoidable and neglect- 
ing thermal fluctuations, Fk remains finite in the limit of 
zero vq. Instabilities and small- velocity kinetic friction 
in more complex model systems have been investigated 
since the Prandtl Tomlinson model was introduced. In 
particular, the Fk(t>o) relationship of elastic manifolds 
sliding in ordered and disordered media has been stud- 



ied extensively 
monly found to be 



The functional dependence is com- 



F k (v)=F h (0)+cv fi . 
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While the elastic instabilities considered in those mod- 
els are important in various contexts, this does not seem 
to be the case for the atomistic explanation of solid fric- 
tion: Many detailed calculations and atomistic computer 
simulations reveal that in almost all cases, interbulk in- 
teractions are too weak to lead to instabilities at the 
atomic scale and as a consequence -Ffc(O) = ||. If inter- 
bulk interactions are very strong, irreversible processes 
like plastic deformation, material mixing, cold-welding, 
etc. usually occur and prevent the instabilities from be- 
ing elastic. 

It has been suggested that the presence of adsorbed 
particles, i.e. a boundary lubricant, confined between 
two surfaces is a more likely explanation for the com- 
monly observed presence of solid friction ^|. The main 
argument is that molecules that are only weakly bound 
to either surface can accommodate the surface corruga- 
tion of both walls simultaneously, which locks the walls 
together. This argument leads to static friction, which is 
the minimum force to initiate sliding between two solids. 
However, in order to explain kinetic friction, it is neces- 
sary to analyze the adiabatic solution x a< \ (i) of the lubri- 
cant atoms. Xgd(t) is the athermal trajectory of an atom 
that always relaxes to the closest mechanical equilibrium 
position at every instant of time. The analysis of x a d(t) 
has become common practice in the context of the mo- 
tion of elastic manifolds in disordered media || . In the 
impurity limit, interactions between lubricant atoms can 
be neglected and consequently x a( j (t) merely depends on 
the initial condition and the relative motion of top and 
bottom wall. In the following, x a d(i) and its connec- 
tion to kinetic friction will be examined. To the best 
of the author's knowledge such an analysis has not yet 
been done, although the present model has already been 
used extensively to predict successfully various tribolog- 
ical phenomena p0 1 1 1 . 



1 



The equation of motion for a lubricant atom in the 
boundary regime reads: 

x = -j b x - 7 t (i - v ) + —T(t) 

m 

--^-{V b (x) + V t (x~v t)}, (2) 

771 OX 

where x denotes the atom's position, m is the atom's 
mass, 74 and 7& parametrize the damping forces from the 
top and the bottom wall, and V t and V& denote the inter- 
action of the confined atom with slider and the substrate. 
r(t) is a Langevin type stochastic random force defining 
temperature. For the lubricant wall interactions, various 
choices will be considered. Centrosymmetric choices can 
all be written in the form: 

Vt. b = V t f cos(avV) + V t f cos(2x/6 t , b ) + ... (3) 

where 27r6 t and 2nbb are the period of the top and 
the bottom wall, respectively. The relevant physical 
units are defined through the choice v} ^ = 1, m = 1, 
bt = 1, and Boltzmann's constant ks = 1- Further- 
more, we will only consider slightly underdamped dy- 
namics (7b = 1 => x(t) w x a d(t)) and restrict ourselves 
to the symmetric choice of Vq := vft , . The free param- 
eters are thus the degree of lattice mismatch (b t — 6&), 
the value of the first higher harmonic (Vi = VSV), and 
temperature T. 

We start the discussion of a commensurate (com.) 
system (b t — &&) without higher harmonics in 
the absence of thermal fluctuations. In that case, 
the net time-dependent potential is simply given by 
V{t) = 2V cos (±vt) cos [l(x - vt/2)] . Thus, for times 
cos (-^vt) ^ 0, the atoms move at velocity vq/2 as shown 
in Fig. [ll An infinitely small moment after a situation 
where this inequality was not satisfied, an atom will not 
be able to find a stable position in the immediate vicinity 
of the previous stable position. As the slider moves on 
with respect to the substrate, the atoms should there- 
fore slide rapidly towards a new mechanical equilibrium. 
However, in order for this to happen, one needs a sym- 
metry breaking element such as thermal fluctuations, 
76 7^ 7t) bt 7^ bb, or round-off errors in numerical cal- 
culations. Otherwise all stable trajectories simply are 
x = titt + 7j t/2 where n is an integer number and the 
average (conservative) force of the lubricant atom act- 
ing on either wall will average to zero. For the athermal 
com. system, we chose to break symmetry by violating 
Galilei invariance and set 74 = 0. This choice is arbi- 
trary but convenient, and is a better-controlled proce- 
dure than relying on round-off erros. We note also that 
for inc. walls, there is no physical reason for assuming 
the equality 74 — jb- 

Owing to broken symmetry (74 = 0), the atom can now 
slide to the next minimum indicated by the gray lines in 



Fig. |l|b. However, the peak velocities x p in this pro- 
cess remain much smaller than during a pop within the 
Prandtl Tomlinson model. The reason is that due to the 
symmetry of Vt and Vb, the pops occur between equiva- 
lent positions in the limit of arbitrarily small v$. Since 
no lower bound for the dissipated energy can be given, 
the zero- velocity Fk will be zero. At the same time, there 
is no upper bound for the ratio x p /vq, so that Fk cannot 
simply vanish linearly with vq but only with some power 
(3 smaller than unity. 
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FIG. 1. Mechanical equilibrium positions for adsorbed 
atoms between two commensurate solid surfaces as a func- 
tion of the relative displacement Aa; wa n between the walls. 
The grey lines indicates the solutions of Eqs. (^),@ if the 
walls are in slow relative sliding motion, (a) First higher 
harmonic Vi < 0. (b) Vi = 0. (c) Vi > 0. 

The nature of the instabilities changes qualitatively 
when the first higher harmonic is different from zero. 
Hence, in the sense of Morse theory |l2|], which contains 
Landau's theory of phase transition as a special case, the 
com. system without higher harmonics can be considered 
a (multi-/tri-) critical point. If V\ < 0, the adiabatic 
solution x(t) becomes continuous as shown in Fig. [I], 
however, the time derivative x(t) diverges. If V\ > 0, 
the adiabatic solution is discontinuous and the pops oc- 
cur between inequivalent positions as shown in Fig. [l]c. 
In analogy to phase transitions, the pops occurring for 
Vi > shall be called first-order instabilities, those for 
Vi < second-order instabilities. Only for first-order 
instabilities can one expect finite energy dissipation and 
thus finite Fk as v$ approaches zero. A numerical analysis 
shows that all three cases can be described with Eq. (Q) . 
The results are shown in Fig. || and the expected trend 
is confirmed: The more discontinuous the adiabatic so- 
lution, the larger the kinetic friction. It is important to 
emphasize that the exponent (3 (as determined for suffi- 
ciently small tjq) only depends on the sign of V\ but not 
on its precise value as long as | V\ \ is not too large. 

For 1-dimcnsional, inc. surfaces, the basic picture is 
similar. If V\ is larger than a (positive) critical value 
Vi , whose precise value depends on the lattice mismatch, 
then pops between inequivalent positions are present and 
Fk remains finite as vo tends to zero. For Vi < V* , 
however, the time derivative of the (continuous) adia- 
batic solution remains finite at all times. Hence the 
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microscopic (peak) velocities v p scale linearly with Vq, 
which implies Stokes-type friction in that regime (1-d, 
inc.). Again, the exponent {3 only depends on the sign of 
Vi — V{ . More detailed results including an analysis of 
significantly more comples systems will be presented in a 
separate paper 
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FIG. 2. Kinetic friction force Ft as a function of sliding 
velocity Vq for commensurate walls with different first higher 
harmonics. Straight lines are fits to low vq data according to 
Eq. (|l]). The two data sets with Vi < are fitted with the 
same exponent /3 = 0.74(0). 

There is however an important difference between com. 
and inc. systems, when the analysis is extended to two- 
dimensional interfaces: In 2-d inc. systems, first-order 
instabilities can now occur without higher harmonics. 
The reason is that in 2-d, atoms can circumnavigate the 
points of maximum longitudinal force. In 2-d com. sur- 
faces, the behavior remains qualitatively similar as in 1- 
d, because the interference of Vt and VJ, remains similar. 
Large-scale molecular dynamics simulations by He and 
Robbins support this argument even if the boundary lu- 
bricant is not any longer in the impurity regime. They 
find small Fk and large F s between com. surfaces, while 
no such gap is observed for inc. systems @]. Also exper- 
imentally, signs for the effects of increased static friction 
of com. surfaces in the presence of a lubricant were re- 
ported [p|. 

This has potentially measurable implications for the 
transition from stick-slip motion to smooth sliding. In 
the stick-slip regime, friction is dominated by F s , while in 
the smooth sliding regime, only Fk is relevant. We have 
extended previous simulations to support this point 
further. For this purpose the same model has been em- 
ployed as that in Ref. || . A schematic of the simulation 
is shown in Fig. || together with the average friction, as 
defined by the energy dissipated per slid distance. The 
verification of this prediction requires smooth surfaces, 
because rough surfaces automatically lead to inhomoge- 
neous energy landscapes The data shown in Fig. || 
was produced at a normal pressure of 0.4 GPa and a ve- 
locity of about 1 m/s using the same conversion of units 
as in Ref. M . Note that the exceedingly small inertia of 



the slider as compared to experiment moves the transi- 
tion from stick-slip to smooth sliding to large velocities. 

We will now turn to the discussion of the effects of 
thermal fluctuations. He and Robbins found that veloc- 
ity dependent corrections in a 2-dimensional, lubricated, 
inc. interface satisfy [13 ]: 

F k M = F k (v Ief ) + O(ln(v /v Ief )) (4) 

over several orders of magnitude in vq. This is different 
from the behavior found in the Prandtl Tomlinson model. 
Although similar temperature corrections have been sug- 
gested, for instance by Prandtl himself, more rigorous 
treatments yield corrections of order (Tin?;) 2 / 3 |jl8| . 
While these corrections describe atomic force microscope 
experiments of nanoscale single-asperity contacts fairly 
accurately, it seems that they cannot provide an explana- 
tion of the usually observed In vq corrections in a straight- 
forward manner. 
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FIG. 3. Kinetic friction force Ft devided by load L as a 
function of spring constant k for com. and inc. walls lubri- 
cated by a quarter layer. A schematic of the simulation is 
shown as well. 

Here it will be shown that simple logarithmic correc- 
tions are obeyed even in the impurity limit. Moreover, 
the crossover to linear response at extremely small slid- 
ing velocities will be included in the discussion. Fig. 3 
shows the normalized friction force obtained at thermal 
energy T = 0.07 for the com. V\ > model. Three 
regimes can be identified. At very small velocities, fric- 
tion is linear in vq and one may associate this regime 
with the creep regime. At intermediate Vq, Eq. (^) is 
rather well satisfied. At "large" velocities, thermal fluc- 
tuations become less relevant and the motion is close to 
that of the athermal system. The data obtained at dif- 
ferent temperatures can be collapsed on a master curve. 
The collapse requires two dimensionless scaling functions 
r(T) and s(T) that both depend on temperature T only. 
The collapse is done via F k (v, T) = r(T*)F k (v*,T*)/r(T 
with s(T)\n(v*(T*)/vi) = a(T*)\n(v(T)/vi) where v x 
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is a constant, r(T) is almost constant [r(0.02) = 
0.98, r(0.2) = 1.2] and s(T) w k b T/AE where AE can be 
interpreted as an effective (free) energy barrier. Qualita- 
tively similar crossover from the linear response regime to 
the activated regime are observed in many other systems 
such as single particles in a static periodic potential [jl9| , 
driven thermal elastic manifolds |^0| , and shear-thinning 
fluids 0. 



1.0 






0.8 


• T = 0.03 
o T = 0.07 
xT = 0.12 


athermal J? 
regime § 


0.6 


vT = 0.2 




0.4 
0.2 


linear response 
in thermal jjP 
- equilibrium,/^ 2 ^ 




0.0 


_K _CL _A 


_Q _Q _1 



10" 6 10" 5 10" 4 10~ 3 10" 2 10" 1 10' 

V 



FIG. 4. Scaling plot of kinetic friction Fj*(v) normalized 
by athermal zero- velocity limit F s — Fk(v = 0, T — 0). The 
reference temperature in this plot is T — 0.07. In the regime 
of thermal equilibrium a linear law Fk(v) oc v is drawn to 
guide the eye. The critical v" contribution is subtracted from 
all data. 

In conclusion, this Letter provides a classification 
scheme for instabilities that can occur when two solids 
are in relative sliding motion. First-order instabilities are 
defined as pops of atoms (or other degrees of freedom) 
between inequivalent positions. They lead to kinetic fric- 
tion that remains finite when the sliding velocity vq goes 
to zero provided the system is athermal. The exponent 
f3 of the velocity corrections v" depends on the details of 
the model, however, p only changes its value at certain 
critical points in the parameter space defining interac- 
tions and geometry. At finite temperature, corrections 
in the order of In vq apply. Second-order instabilities are 
defined as pops between equivalent positions. They lead 
to a sub- linear power law oc v" . If, however, the adia- 
batic solution x a d of the boundary lubricant moves with 
finite velocity at all times, then simple Stokes friction 
follows. Thus an important result of this analysis is the 
identification of critical points where the kinetic friction 
law changes qualitatively as a parameter describing the 
interactions or the geometry is varied. 

While the present study is primarily concerned with di- 
lute boundary lubricants, the concept itself seems to be 
more general. For instance x a d can be a collective order 
parameter that fluctuates back and forth between two 
values. Such quasi-periodic phase transitions have been 
observed in computer simulations of Ni asperities moving 
over a Cu substrate [E2J. Of course, the situation in those 



simulations was more complex, because wear occurred as 
a side effect of the motion. In the other extreme, a; a d 
may merely denote the position of an electronic orbital. 
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